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The coherent control of scattering processes is considered, with electron impact dissociation of 
used as an example. The physical mechanism underlying coherently controlled stationary state 
scattering is exposed by analyzing a control scenario that relies on previously established entangle- 
ment requirements between the scattering partners. Specifically, initial state entanglement assures 
that all collisions in the scattering volume yield the desirable scattering configuration. Scattering 
is controlled by preparing the particular internal state wave function that leads to the favored col- 
lisional configuration in the collision volume. This insight allows coherent control to be extended 
to the case of time-dependent scattering. Specifically, we identify reactive scattering scenarios us- 
ing incident wave packets of translational motion where coherent control is operational and initial 
state entanglement is unnecessary. Both the stationary and time- dependent scenarios incorporate 
extended coherence features, making them physically distinct. From a theoretical point of view, 
this work represents a large step forward in the qualitative understanding of coherently controlled 
reactive scattering. From an experimental viewpoint, it offers an alternative to entanglement-based 
control schemes. However, both methods present significant challenges to existing experimental 
technologies. 



I. INTRODUCTION 

Coherent control [l[ is an approach to controlling quan- 
tum processes where the phase coherence between dif- 
ferent quantum states is explicitly used in order to en- 
hance or suppress a desired outcome of a given quantum 
event through interference effects. In the case of scat- 
tering, control over the product cross sections has been 
shown to be achievable by creating a coherent superpo- 
sition of incident scattering states @, [|[ . Unlike the co- 
herent control of unimolecular processes, which has been 
established both theoretically and experimentally |l[ , ap- 
plications of coherent control to collisions is only in its 
infancy. Hence it is important, at this early stage, to 
clarify and reinforce its essential principles, as well as to 
identify and develop qualitative pictures and insights ap- 
plicable to a wide class of scattering problems, including 
the all-important case of reactive scattering. 

Previous studies focusing on crossed beam scenar- 
ios at fixed total energy have identified necessary require- 
ments for the coherent control of scattering events. In 
these stationary state scenarios, incident states consisting 
of coherent superpositions of internal and translational 
motions were considered. In general, due to the conser- 
vation of center-of-mass momentum and of energy during 
the collision, initial state entanglement between the inter- 
nal and center-of-mass states of the incident beams was 
found to be required. Here entanglement refers to the 
non-separability of the initial wave function when writ- 
ten in terms of the lab frame coordinates of the incident 
scattering particles. This requirement presents a consid- 
erable challenge to experimental implementation of the 
coherent control of scattering processes, and is partially 
responsible for the lack of coherent control scattering 



experiments thus far. The need for initial state entan- 
glement in stationary state scenarios can be removed if 
superpositions of degenerate internal states are usedQ. 
However, reliance upon such entanglement-exempt sys- 
tems greatly restricts the choice of possible situations 
that can be used in the coherent control scenario. 

Although the formalism leading to conditions for the 
coherent control of stationary state scattering is clear @, 
the associated qualitative insight into the role of initial 
state entanglement is lacking. In this paper, we expose 
the universal configuration-based mechanism that lies at 
the core of such coherently controlled scattering Al- 
though our results are applicable to all scattering pro- 
cesses, we focus below on the most challenging case of re- 
active scattering. In particular, by analyzing the generic 
fixed total energy control scenarios where entangled in- 
terfering pathways are required, it is found that entan- 
glement between the two incident beams ensures that all 
collisions within the scattering volume occur at a single 
favorable configuration of internal states. Further, con- 
trolling the relative phases of the entangled interfering 
pathways is shown to shape the internal state wave func- 
tion that participates in the scattering event. 

Having identified this general mechanism, we extend 
this thinking to the time-dependent scattering regime, 
allowing us to introduce an alternate approach to con- 
trolled reactive scattering that uses non-entangled wave 
packets of translational motion. This approach relies on 
the observation that translational wave packets have a fi- 
nite size, so that there is a finite volume of overlap in both 
space and time, defined by the scattering beams, wherein 
the collisions occur. If the internal configuration does 
not have time to change as the molecules move across 
the collision region, then a desired configuration can be 



2 



established and maintained during the collision time, al- 
lowing coherent control of reactive scattering despite the 
absence of initial state entanglement. This method of 
control relies on the presence of temporal correlations 
in the incident state (the collisions need to be precisely 
timed relative to the internal motion) as opposed to ini- 
tial state entanglement. 

As a result of this work we have exposed the underly- 
ing qualitative mechanism for coherently controlled scat- 
tering and developed it within both stationary and non- 
stationary frameworks. Each approach has its own chal- 
lenges in experimental implementation, and each has its 
individual benefits. For example, the time-dependent 
wave packet version is perhaps conceptually and intu- 
itively simpler, while the stationary state entangled ver- 
sion is not limited to scenarios with restricted collision 
volumes in space and/or time. 

In the molecular scattering literature, it is frequently 
the case that one uses time-dependent (i.e. wave-packet) 
methods to calculate essentially time-independent scat- 
tering properties [H, such as energy-resolved cross sec- 
tions. That usage of wave packets is intended solely 
as a computational convenience, and the same scatter- 
ing results could have been computed using fully time- 
independent methods. We emphasize that, by contrast, 
the stationary and time-dependent scattering control 
schemes that we discuss do not simply refer to two dif- 
ferent computational methodologies used to calculate the 
same physical scattering scenario. Rather, the stationary 
and time-dependent control scenarios embody physically 
distinct extended coherence properties: one requires ini- 
tial state entanglement while the other requires temporal 
correlations between the incident wave packet. 

Similarly, the time dependent scheme introduced be- 
low does not just consist of simple shaping of a transla- 
tional wavefunction incident on a single state of a target 
molecule. Such a scenario would be akin to the previ- 
ously studied 2PACC scenario Q which we have recently 
shown[6[ does not involve any aspect of quantum inter- 
ference, and is hence not coherently controlled collision 
dynamics. Rather, as shown below, our time dependent 
approach invokes specific interference effects to affect the 
product cross sections. 

As a working example to illustrate these ideas, we 
consider electron impact dissociation of Hj, where the 
total dissociation cross section and energy spectrum of 
the ejected protons are the observables to be controlled. 
However, it should be stressed that the arguments and 
conclusions to follow are completely general and could 
have been illustrated using any scattering example. Sim- 
ilarly, none of our results relies on the use of the Born 
approximation that is invoked below. However, the spe- 
cific example of e-Hj collisions is focused upon since it 
is of central interest to the emerging field of strong field 
attosecond physics Q , where laser-induced ionization fol- 
lowed by electron recollision with the parent ion lie at the 
core of many processes. 

The paper is organized as follows: Section HT1 provides 



a reformulation of controlled stationary state scattering 
which affords insight into the origin of conditions for con- 
trol, and Sect. |TTT] displays the underlying mechanism for 
control via computations on electron impact dissociation 
of H j" . The extension of these qualitative control insights 
to the case of time-dependent scattering is discussed in 
Sect. IhTbI Section ITVl provides a summary. 



II. FORMALISM 

A. General stationary state scattering 
considerations 

Typical crossed molecular beam experiments are well 
described by stationary state scattering theory at fixed 
total energy. Consider then a few general results from 
stationary state scattering theory. Within the S'-matrix 
formalism [8, 9] the transition probability from the initial 
la) to final lb) state is 
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where H is the scattering Hamiltonian, T a b is the transi- 
tion matrix, and (5 a b represents the unscattered compo- 
nent. Note that all equations are written in atomic units 
h = m P = e = 1. The T- matrix elements have the form 



T ab = 2ir 5{E h - E a ) <5(K a - K b ) t ab , 



(2) 



where E a and E^ are the energies of the initial and final 
states, K a and Kb are the center-of-mass momentum of 
these states. The i a b are the "on-shell" T-matrix ele- 
ments, and depend only on the relative momentum, k a 
and kb, and the internal quantum numbers, denoted a' 
and b', of the incident and outgoing states. Physically, 
the two (5-functions in Eq. ((2|) enforce conservation of 
energy and of total center-of-mass momentum. 

When the scattering pair is launched in a single eigen- 
state | a) of the reactant system, the differential cross 
section is 



d<7b' (a) 
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(2n) 4 (l/|k a |)f 
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(3) 



where Kb = K a , the magnitude of the outgoing relative 
momentum is 



|k b | = V2(|k a | 2 /2 + £a< -E h ,), 



(4) 



and SI is the 3D angle of the vector kb- The total cross 
section for a particular reactive arrangement channel n is 
then found by summing over all internal states belonging 
to n, and integrating over the solid angle O, 
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Since the on-shell transition matrix elements natu- 
rally depend on the initial state |a), one way to en- 
hance/suppress a desired reactive cross section is to use 
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a single incident state |a), and vary it until the particu- 
lar incident state is found that achieves this goal. This 
method is termed passive single-state control. 

Experimentally creating a single incident scattering 
eigenstate is not always an easy task. Often, internal 
degrees of freedom of the scattering pairs are in a ther- 
mal distribution of quantum states. In this situation, 
cross sections can be controlled to some degree by vary- 
ing, for example, the internal temperature of the particles 
leading to a temperature-dependent cross section 

^ ( " ) (T) = i^e-^/ fc V")(a'), (6) 

where Z is the partition function. Varying the temper- 
ature then changes the distribution of states that par- 
ticipate in the collision, and hence offers control. The 
thermal case can, of course, be identified as a particular 
case of scattering from an incoherent initial distribution, 
which in general leads to cross sections of the form 
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(7) 



where F(&) defines the incoherent distribution of incident 
states. Since all incoherent distributions provide averages 
over the single-state cross-sections, one can never achieve 
greater controllability using incoherent distributions than 
that achievable in single-state scattering. However, the 
use of coherent superpositions of incident eigenstates, in- 
stead of incoherent distributions, provides an opportu- 
nity to do better than the single-state scenario. 



B. Superposition states and coherent control 



The results in Sect. Ill Al may be extended by consider- 
ing scattering from an initial superposition state, either 
of fixed energy (and hence stationary) or of varying en- 
ergy content (and hence non-stationary). The full differ- 
ential cross section for scattering from the coherent initial 
superposition 



(O. However, full dimensionality is allowed in the outgo- 
ing states. Allowing for off-axis incident momenta would 
simply introduce addition integrals over the initial mo- 
menta to Eq. @ , but the results of the present study are 
otherwise unaffected. The total cross section for scatter- 
ing into channel n is obtained from Eq. 



as 



» 



dkhdKb 



(11) 



Within the coherent control approach, the relative 
phases between multiple pathways from the initial to final 
state are used to control the process. By controlling the 
relative phase of the (complex-valued) C(a', fc a , K a ) in 
the initial state Eq. ([8]), multiple coherent pathways can 
be created that manifest themselves as the sum over a' in 
Eq. From the form of the T-matrix [Eq. ©J, we see 
that all allowed transitions must conserve center-of-mass 
momentum and total energy. This leads immediately to 
the general requirement for scattering interference: the 
initial states contain on-shell coherence, that is, in or- 
der for two initial eigenstates to interfere, they must be- 
long to a single shell as defined by the 5-functions in the 
T a b-matrix. Hence, pathways in Eq.® exhibit on-shell 
coherence (since K a = Kb, and Eq. (TIU|) is simply a 
statement of energy conservation) if the initial superpo- 
sition Eq. ([5]) contains more than one state satisfying 
these on-shell requirements. 

In the case of field-free bimolecular scattering governed 
by Eq. as will be evidenced below, superposition 

states with coherence in the lab frame translational mo- 
tion of both incident particles, and at least one internal 
mode (minimum 3 degrees of freedom in total), is re- 
quired for interference. If static or time-dependent ex- 
ternal fields arc present during the collision event then 
they will modify the on-shell conditions, through energy 
and momentum exchange with the particles, and can lead 
to less restrictive conditions on the required initial-state 
coherence. 



C. e-H^ scattering 
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tron impact dissociation of Hj ■ The Hamiltonian for this 
scattering problem is given by 



dkbdKb 
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(9) 
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where again K a = Kb, and the initial relative momenta 
corresponding to each kb and b' is now given by 



|Ri-? b | |R 2 -r b | 
1 1 



k a = ^/2(\kb\ 2 /2 + E h >-E a 



(10) 



R2 — r Ri — R2 



|Ri - r| 
1 

|r b - r I ' 



For simplicity, the assumption that all incident momenta 
lie along a single axis (i.e. the initial momenta of the 
two particle are antiparallel) was used to arrive at Eq. 



where the momentum/position operators with the sub- 
scripts T' and '2' refer to the two protons, those with 
the 'b' subscript refer to the bound electron, and the 
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remain operators refer to the incident electron. In this 
process, an electron with momentum is incident on an 
molecule of momentum P^, which is in an internal vi- 
brational and rotational state labeled by n = (y, J, mj) 
with energy E n on the ground electronic state £ s . All 
indicated momenta are in the laboratory frame. Dur- 
ing the e-Hj collision, the incident electron excites the 
bound electron from the bonding to antibonding state, 
£ g — > £ u , through the electron-electron Coulomb inter- 
action V ee - The final state of the scattered particles con- 
sists of the scattered electron with momentum pf and 
two protons with momentum Pi and P2, one of which 
carries the bound electron. For the purposes of the scat- 
tering calculation, the final state of the two protons is 
written in terms of the center-of-mass motion of the full 
composite Pf = Pi + P2 and the relative motion 
P r = (Pi — P2V2 corresponding to a continuum state 
of energy E = |P r | 2 /2/i and angular momentum state L 
on the surface of H^, where fi = m p /2 is the reduced 
mass of the molecular ion and m p is the mass of the pro- 
ton. The E u continuum states asymptotically approach 
the free particle momentum states defined by Pi and P2 
at large distances, but are distorted near the core. 

The on-shell transition matrix elements can be eval- 
uated to first order in the electron-electron interaction 
(first Born approximation) [IcL [TlT | , 



(13) 



where the energies are given by 



E b 



2 

r> 2 



2rri] 
2mj 



(14a) 
(14b) 



Kj and kj (j = i, /) are the initial and final center-of- 
mass and relative momenta of the e-H^ system 



K, 



mi + 1 



(15a) 
(15b) 



and mj 
electron-electron 



2m p is the mass of the ion. The 



interaction 



matrix 
evaluated 



elements, 
in full- 



dimensionality using: the LCAO approximation for 

the bound electron [ll| as described in Appendix A. 
We avoid additional approximations [l(| [H| by using 
numerical wave functions for the radial molecular 
continuum states, with the E g and E n surfaces taken 
fromRcf. [l|. 

The full differential dissociative cross section is 



^i D) (*o) 
dEdkfdK 



(27T) 4 



n 

1 



(16) 



where 



III. COHERENT CONTROL OF REACTIVE 
SCATTERING 



k M = J2(\k f \y2 + E-E n ), 



(17) 



and the initial momenta are again assumed to be an- 
tiparallel. The cross section da^ D ' /dE, dependent only 
on the energy of the ejected protons E, is calculated by 
integrating over the unobserved coordinates 



dE 



E 

L=l,3,5... 



dEdkfdK 



dkfdK 



and the total yield is given by 



(18) 



A. Stationary State Scattering: Few-state 
entangled superpositions 



From the requirement of on-shell coherence, the sim- 
plest (in terms of number of states involved) incident su- 
perposition [Eq. ©] that offers on-shell coherence, and 
hence coherent control, utilizes two different a' states, 
both with the same total energy (i.e. internal plus trans- 
lational energy), at the same center-of-mass momentum 
K, i.e. a state of the form: 



r (D) (* ) 



dE 



dE. 



(19) 



Ci|ai)|*i)|Jf) +e^C 2 |a 2 )|fc 2 )|if), (20) 
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where C\ and C 2 are real. The associated differential 
cross section is then 



dh'dK 



(2 7 r) 4 |C 1 (l/fc 1 )t aib , +e^C 2 (l/fc 2 )^ b ,| 2 (21) 



(2tt) 



\Ci(l/h)t a{b >\ 2 + \C 2 (l/k 2 )t^ b > 



+ \C 1 (l/k 1 )t a , ih ,\\C 2 (l/k 2 )t^ b ,\cos(<p < y 2 + <j> 

where W^a!, is the phase of i*/ b /i a ;,b'- Varying the rela- 
tive phase 4> of the two components in Eq. (f2T)|) controls 
interferences in the scattered state and provides a means 
of controlling the scattering products. By tuning <j) one 
can enhance or suppress the cross sections beyond the 
values attainable with single-state and incoherent sce- 
narios that utilize the same states, provided that t a ' ib i 
and t a ^b' are not too dissimilar. 

The state \^ s) is expressed above in center-of-mass 
and relative coordinates. In lab frame coordinates, ex- 
plicitly for the example of e-H^" scattering, \^ os ) becomes 

|*„.) = dlnW)!^)^) + C 2 e i V 2) >l^ (2) >bf ) > ! 

(22) 

subject to the constraints 



2mj 



(1)^2 



d(2) , (2) 
= Pi +Pi 



(23a) 



(2)n2 



(PTT , (pr) 

2m/ 



Note that, in the general case when the internal state 
In^ 1 )) and |n' 2 ') are not degenerate, Eq. (|2"2"|) is an en- 
tangled state of translational and internal motion, as pre- 
viously discussed [H- For degenerate internal states, Eq. 
(gH permit the solution pf } = Pp\ pf ] = pf\ thus 
removing the initial state entanglement requirement Q . 

In the remainder of this section we analyze this con- 
trol scenario to expose the qualitative mechanism lying 
at the core of the entangled stationary state control, and 
then use this insight to introduce a time-dependent non- 
entangled version of coherently controlled reactive scat- 
tering, described further in the following subsection. 

An example of control using the state Eq. (j22f is 
shown in Fig. [TJ Specifically, this example superposes 
two molecular vibrational states with v = and 1 and 
with energies Eq = -0.0973 au and E\ = -0.0871 au, angu- 
lar momentum J = mj — 0, and translational momenta 
pj 1 ' = au and pf' = 4 au The momenta P^ and p\ 2 ^ 
are then set in accordance with Eq. (|23[) . The weights 
of the two components are set equal, C\ = C 2 . Panel 
(a) shows the total cross section for dissociation as <j> is 
varied. For comparison, the two dashed lines show the 
cross sections for the v = (lower line) and v = 1 (upper 
line) states, considered separately. Enhancement or sup- 
pression of the total cross section beyond the incoherent 
result is clearly evident in Fig. [1] when using the entan- 
gled coherent superposition Eq. (|22p . The remaining two 





FIG. 1: Control of reactive scattering using the entangled 
initial state Eq. (|22|l . Panel (a) shows the dependence of 
the dissociation cross section a' 13 ' on <f). The dashed lines 
denote cr' ' for a single component wave function using v = 
(lower line) and v = 1 (upper line). Panel (b) plots the 
energy spectrum of the proton fragments using v — (solid) 
and v = 1 (dashed). Panel (c) shows the analogous spectra 
when using the two-state superposition with (f> — (solid) and 
4> = 7r (dashed). 



panels show the energy resolved cross sections, da^ /dE, 
for the states v = and v = 1 individually [panel (a)], 
and for the superposition states corresponding to the ex- 
trema of (</>), namely <f> = and n [panel (b)]. A 
large degree of control over the proton energy spectrum 
is evident. 

An important physically transparent qualitative pic- 
ture of how the control arises for the entangled state can 
be constructed. To do so, we first note that there is am- 
ple evidence 0, [T§] that the exchange kernels between 
the reactant and product arrangements are short range. 
As such, the behavior of the wave function when the re- 
actants are close to one another is particularly relevant 
to the reactive cross sections. For this reason, we fo- 
cus below on the character of the wave function at short 
range. 

In relative and center-of-mass momenta, the incident 
state is 



C^)\k^)+C 2 e^)\kf) 



\K- 



(24) 



Since the center-of-mass momentum is conserved dur- 
ing the scattering, only the terms in the square brack- 
ets in Eq. (|24| are relevant for the control dynam- 
ics. Figures [2ji and [2)3 show the probability P(R, x) = 
\(R, x\i$> os )\ 2 , of finding \^ os ) at the nuclear bond length 
R and electron-ion separation x (the conjugate of k{). 
The internal states and incident momenta are the same 
as those used in Fig. [TJ panel (a) shows results for = 0, 
while panel (b) uses cj) = ir. Note that these plots show 
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the incident states in the absence of interparticle interac- 
tions; they do not include the scattered component. The 
addition of the latter will change our argument quantita- 
tively but not qualitatively. As noted above, of particular 
interest is the character of the internal state near the col- 
lision region, x ~ 0. Since \^ff s) is a time-independent 
wave function the probabilities plotted in Figs. [5^, and 
reflect the complete incident dynamics of the scat- 
tering pair. Panels (a) and (b) show that the internal 
vibrational wave function near x = 0, in fact for all x, is 
controlled by <j>. Further, from Eq. (pM|) . it is seen that 
the wave function at x = is 



(R,x = Q\V os ) = \aiR\vW) +C i e*(R\v&>) 



(25) 



which, by varying <f>, can be shaped into structures not 
accessible using individual incident eigenstates. Since the 
cross section is strongly dependent on the internal wave 
function, controlling (j) then allows one to control en- ) 
and dcr( D '/dE by manipulating the particular internal 
configuration that participates in the collision. This is 
the qualitative mechanism lying at the core of coherently 
controlled reactive scattering using the entangled super- 
position state Eq. (f2"2"|) : The initial superposition state 
defines the structure of the internal state wave function 
at the point of collision and can therefore be used to tune 
this structure in order to optimize the reactive cross sec- 
tion. 

In standard scattering theory one often considers 
scattering of incident time-independent eigenstates, and 
hence this observation may seem trivial. This is not the 
case. Had a non-entangled wave function of the form 

l^ir') = [C 1 |n( 1 ))+C 2e ^|n( 2 ')] (26) 
x [C 3 |^ (1) )+C 4 e*|^ (2) )] 
x [C 5 \p^)+C 6 e^\p^)} 



been used, the resultant incident state would contain 8 
terms, only two of which are at the same total energy 
and center-of-mass momentum, namely the components 
in Eq. (|22"|) . From an energy domain perspective, this 
means that only two of the total eight terms exhibit on- 
shell coherence, and hence only two of the eight terms 
can interfere. The remaining six terms contribute in- 
coherently to the cross sections. From a time domain 
perspective, only the superposition of the two on-shell 
components lead to a time-independent wave function. 
All other components, since they are at different ener- 
gies, accumulate a time-dependent phase relative to the 
on-shell superposition. With respect to the internal con- 
figuration near the collision region, this means that the 
off-shell components will alternate between constructive 
and destructive interferences at different times, introduc- 
ing a time-average over all possible phases relative to the 
on-shell component, and hence an average over internal 
configurations allowed by the populated internal states at 
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FIG. 2: Plots of the probability density of relative and internal 
motion P(R,x) = \{R,x\^)\ 2 . initial state Eq. <(22j). Panels 
(a) and (b) correspond the case shown in Fig. Q] for <j> — 
and n respectively. Panels (c) and (d) show incident states 
constructed by hand to minimize and maximize the nuclear 
bond length at the moment of collision, and hence minimize 
and maximize a^ D \ 



the collision point. These time- varying interferences re- 
duce the off-shell contribution to an effective incoherent 
contribution. The entanglement in the few-state super- 
position Eq. (|22[) then plays a crucial role in limiting the 
particular internal state wave functions that participate 
in the scattering events. 

Further investigation of the e-H^" cross section 
strengthens the configuration-based qualitative mecha- 
nism described above. In particular, it is known that 
the electron impact dissociation cross section of is a 
strongly dependent monotonic function of bond length i?; 
as (R) increases, so does [14| . Comparing the x = 
structure in Figs, [^i and [2Jd shows that the control of 
<j( D > follows the expected (-R)-dependence; the <j) = tt 
case has both the larger (R) and a^ D \ 

Using this insight, a clear extension of the two-state 
entangled coherent control scenario to a multi-state en- 
tangled scenario, in order to achieve better enhance- 
ment/suppression than that of Fig. [TJ can be con- 
structed. Since control here is directly linked to the aver- 
age bond length at x — 0, the superpositions that maxi- 
mize/minimize the scattering cross sections are the same 
superposition that maximize/minimize (R) at x — 0. For 
example, two superpositions (constructed by hand) that 
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approximately achieve these goals are 



l^max) — -^max 



and 



|*min) = N a 



j2(-iye-(^y\v)\k v ) 



\K) (27) 



E 1 



\K) (28) 



where iV max and iV m i n are normalization constants, fco 
= 4 au, and the remaining fc„ are set to ensure that 
all components satisfy the on-shcll requirement, hav- 
ing the same total energy and center-of-mass momen- 
tum. The summation runs over all the vibrational states 
The corresponding P(R,x) are shown in Figs. &p 
and [2ji. For these states, the values of (R) in the 
collision region are ( x I > max|-R|^'max)x=o = 8.75 au and 
(^'min|-R| x I > min)x=o = 1-32 au, while the largest and small- 
est values possible using an incoherent initial state would 
correspond to those of the highest and lowest vibrational 
states of H2 , i.e., {v = lS\R\v = 18)2=0 = 8.12 au and 
{v = 0|i?|z/ = 0)a,. = o = 2.05 au The coherent superposi- 
tions l^max) and I'i'min) access values of (R) at x — be- 
yond those accessible to any incoherent mixture of states, 
and hence provide more control over than any inco- 
herent scenario. 

Note that although in the case of e-H^~ scattering there 
is a clear classical explanation of the nature of the con- 
figuration that enhances control, this is need not be the 
case in general. Rather, given the Hamiltonian H and the 
associated scattering matrix S^b, there is a well defined 
scattering configuration that maximizes control. When 
the on-shell (and hence entanglement) requirements are 
met for the initial superposition [Eq. ([5])] then the vary- 
ing the coefficients C(a', k & , K a ) in Eq. ([8]) alters the sta- 
tionary spatial configuration and hence alters the cross 
sections. The optimal choice of the C(a', /c a , K & ) then 
corresponds to the superposition that comes closest to 
the optimal stationary state configuration for scattering. 
Whether the optimal configuration is easily understood 
classically, however, depends upon the system under con- 
sideration. 



B. Non-entangled wave-packet superpositions: 
Time dependent scattering 

Having exposed the qualitative mechanism that un- 
derlies control resulting from stationary state, and hence 
entangled, initial superposition states, additional ap- 
proaches to coherent control of reactive scattering that 
do not rely on initial state entanglement can be identi- 
fied. For example, there is an alternative route to en- 
suring that collisions between two particles occur at a 
fixed phase of the internal-state motion: design time- 
dependent superposition states (i.e. Eq. ([5]) with non- 
degenerate states), specifically wave packets of transla- 



tional motion plus superpositions of internal states, in or- 
der to localize the collision partners in space, and thereby 
to restrict the duration of the collision between the wave 
packet to less than the internal state motion (see Fig. 
[3^,). Quantum interference manifests in this case due to 
the numerous energetically degenerate sets of states that 
occur due to the energy widths of the two incident wave 
packets. 

The initial state used to illustrate the wave packet sce- 
nario is 



where 



\il>u)=[\" = 0)+e i +\v = l)] /yfi, 



dpi(A pi y/Tr) 2 g 



\Pi 



and 



= / dP^Ap^T 



\pi) 



(29) 



(30a) 



(30b) 



(30c) 



as depicted in schematically in Fig. [3^,. Note that this 
is not simply wavepacket scattering off of a single vibra- 
tional state. Rather, scattering is off an internal super- 
position of states, necessary to incorporate interference 
and achieve coherent control. 

To get an idea of how long the collision between the 
two wave packets lasts, the following measure, here called 
the time-dependent collision probability W c (t), is used 



W c (t)=N\(*(jb)\6(x-y)\*(t))\ 



(31) 



where x and y are the positions of the electron and ion re- 
spectively, \^(t)) is the incident wave function (i.e. scat- 
tered components are not included) and N is a normal- 
ization constant such that J W c (t)dt = 1. This quantity 
gives the probability of finding the electron and the ion at 
the same position in space, and hence reflects the prob- 
ability of a collision occurring at time t. In the continu- 
ous beam scenario, as considered in the previous section, 
W c (t) is a time-independent constant. However for wave 
packet collisions, W c (t) will be a localized Gaussian-like 
function indicating that, in this case, collisions only oc- 
cur during a select time window where the two colliding 
wave packets overlap in space. The duration of the wave 
packet collision is then defined as twice the standard de- 
viation of W c (t) 



AW, 



= 2^(*) Wc (t) 



2 



(32) 



where (■ ■ - )w c indicates an average value where W c {t) is 
used as the distribution function. 

In order to compare the wave packet scenario to the 
entangled state scenario, we set the mean incident mo- 
menta of Eq. (|2"2")) to Pio = and pio = 4 au, the same 
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t < t 




t = t 




t > t 





curves. The second method offsets the spatial and tem- 
poral focus of the electron wave packet relative to the ion 
wave packet 



Position 



FIG. 3: Control of reactive scattering using the wave packets 
of translational motion. Panel (a) shows schematically the 
wave packet scenario: controlling the internal state of the 
molecular ion during the finite duration of the collision allows 
for control. t c denotes the time of the collision. Panel (b) 
shows the degree of control using the initial Eq. (|29p . See the 
text for the specific parameters used. Panel (c) demonstrates 
the loss of control as the duration of the collision between 
the two wave packet approaches the timescale of the internal 
state motion. The dashed line denotes the vibrational period 
of H+, T V ib = 14.9 fs. 



as used above. Using then momentum widths of Ap z = 1 
au and A Pi — 0.01 au gives a scattering scenario where 
AW C = 0.87 fs. Since AW C in this case is much smaller 
than the vibrational period of Hj (t v h, = 14.9 fs), we 
expect that control is possible using these parameters. 
Figure [3J3 plots the corresponding reactive cross section 
as 4> is varied. Control is indeed present. Further, upon 
comparing Figs. Wp and [TJi, one sees that the control 
using the non-entangled state Eq. (f29|) gives the same 
degree of control as the entangled state Eq. con- 
firming that the wave packet scenario offers an alternate 
and equivalent (in terms of the controllability of the to- 
tal cross section) means of coherent control of reactive 
scattering, at least for reactive scattering absent of any 
resonances related to the relative momentum fcj. This 
latter point will be revisited in the following section and 
is at the root of possible advantages of using the entan- 
gled scenario. 

The remaining panel of Fig. [3] plots the minimum and 
maximum values of the cross section, which correspond to 
= and 7T, as the duration of the wave packet collision 
AW C is increased. Two different methods of increasing 
AW C were explored. For the first method, Api is gradu- 
ally made smaller and smaller, which causes the spatial 
size of the electron wave packet to increase and thereby 
increases AW C . The control results using this technique, 
plotted against AFT C , are shown in Fig. [3J: using the solid 



\^Pi) = J dpi(A Pi ^) 



(33) 



where Xd = PioTd, and Td controls the focusing offset. 
The control results using this method are plotted in Fig. 
Hfc with circles. In both cases, the control goes to zero as 
the AW C approaches the vibrational period t v h,. These 
calculations most clearly demonstrate the configuration- 
based mechanism underlying the control; the controlla- 
bility decreases as more internal configurations partici- 
pate in the collisional events. 



C. Entangled or Wave Packets? 

Both the stationary state and non-stationary state 
control schemes present technological challenges. The 
wave packet case requires, experimentally, a collision over 
very short times, and therefore must be run in either 
a pulsed mode, or with very tightly focused molecular 
beams with small overlap region such that the molecules 
move through the collision region in a time smaller than 
the timescale characteristic of the dynamics of the inter- 
nal superposition. The stationary state case, on the other 
hand, requires initial state entanglement, but can be used 
in a continuous beam regime and/or with arbitrarily large 
spatial region of overlap of the two beams. Hence, the 
entangled version offers arbitrarily large space-time scat- 
tering volumes while spatial restrictions exist for the wave 
packet version. This implies that, although the cross sec- 
tions may be controlled to the same degree, the entangled 
version will always permit larger total yields since it came 
be used in conjunction with arbitrarily large volumes and 
incident fluxes. 

A second advantage of the entangled scenario relates 
to the possibility of scattering resonances. The reactive 
cross section for e-H^~ scattering is a rather smooth func- 
tion of the incident relative momentum k{. Hence, our 
sample cases thus far have implicitly considered the con- 
trol of reactive scattering in the absence of sharp reso- 
nances related to the incident kinetic energy. However, 
reactive scattering problems of chemical interest may ex- 
hibit such resonances (e.g. Feshbach resonances). The 
pulsed wave packet scenario requires a broad superposi- 
tion of incident momenta in order to obtain spatial local- 
ization and short-time overlaps of the two incident wave 
functions, while the entangled case can use as few as two 
incident momenta. It may very well be the case that nar- 
row resonances can be resolved/exploited in the entan- 
gled case .10] , while averaging over the broad momentum 
bandwidth in the wave packet case washes out the nar- 
row resonance features, rendering them unusable in this 
latter scenario. 
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Note that both of these aspects of initial state entan- 
glement, the possibility to efficiently exploit narrow res- 
onances and access arbitrarily large scattering volumes, 
represent non-classical aspects of controlled reactive scat- 
tering accessible via entanglement. 



(NSERC) of Canada. 



APPENDIX A: ELECTRON IMPACT MATRIX 
ELEMENTS 



IV. SUMMARY 

Previous work on stationary state coherent control, 
which provides a useful description of crossed beam ex- 
periments, was shown to require initial state entangle- 
ment between the incident translational and the inter- 
nal states of the scattering partners @. This paper has 
provided physical insight into the role of this require- 
ment, generating an extension to coherent control via 
time-dependent wave packet scattering. Specifically, the 
initial state entanglement was shown to assure that all 
collisions occur at a fixed configuration of the internal 
state motion. This qualitative insight then allows for 
the introduction of a coherent control scenario in time- 
dependent scattering. Control in the latter regime is pos- 
sible if the duration of the wave packet collision is much 
smaller than the characteristic timescales of motion of 
the superposition of internal states. 

The mechanism of fixed-configuration scattering un- 
derlying coherent control of reactive scattering can be 
extended to all scattering scenarios. A sample extension 
to scattering off surfaces is provided elsewhere @. 

Some extensions of this approach are worth noting. 
First, scattering studies carried out on loosely bound 
van der Waals complexes are distantly related to the 
wave packet control scenario. In these studies, a CO2 • 
HBr complex was used as an oriented precursor to study 
the CO2 + H reaction, where photodissociation of HBr 
launched the H toward the CO2. These experiments are 
examples of the selection of well defined angular wave 
packets of the scattering partners. However, for active 
control of reactive scattering in these systems, one would 
also need to tune the particular angular wave packets 
that participate, as opposed to selecting the single ori- 
entation defined by the initial van der Waals complex. 
This could perhaps be accomplished by selectively excit- 
ing rotational states of the CO2 and/or HBr prior to the 
initiation of the reaction. 

Second, although not explored in this paper, both 
the entangled and wave packet scenarios explored herein 
could be used instead to completely characterize the scat- 
tering matrix 5 a b with phases, in analogy with methods 
of Quantum Process Tomography [l7| . In short, hav- 
ing selective control over the input internal state, in ad- 
dition to the usual control over the input translational 
momenta, would allow one to measure enough projec- 
tions of the scattering operators, through measurements 
of the scattering cross sections for different internal su- 
perpositions, to allow for an accurate reconstruction of 
the complex S-matrix. 
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For both the strong field and field-free scenario, the 
transition matrix elements (^ouil^eel^in) are required. 
In calculating these values, we use the LCAO approxi- 
mation for the bound electron [l(| El ■ Specifically, using 
a basis of definite energy and angular momentum for the 
final state of the internuclear coordinate R and restrict- 
ing the initial Hj to zero angular momentum gives 



% , k, \V ee \ ^ , ki) = i L y/(2L + l)K(L, u,E,k), 



(Al) 



where k = k/ — k;, k = |k|, and 



Tl(L,v, E,Ak) 
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7T 2 k 2 [4 + k 2 } 2 

(R)N^\r) X e,j (R)jl (kR/2) X v{R)dR. (A2) 



The z-axis of the angular states lies along Ak, and only 
the niL — sub-levels along this axis have non-zero am- 
plitude. The N^\R) are normalization factors arising 
from the LCAO wave functions of the bound electron and 
are given by 



N (±) (R) = [2± 2e~ R (l + R + R 2 /i)\ 



-1/2 



(A3) 



The Xv{R) are the bound radial eigenstates on the £ g 
surface 



1 d 2 
^dR 2 



V Sg (R)-E v 



x»{R) = 



(A4) 



while the Xe,l(R) are the continuum radial eigenstates 
on 



1 92 ,y (m 1 L(L + 1) F 



(A5) 

The Vs g (R) and Vs u (R) surfaces are taken from Bunkin 
and Tugov [13[ . Vs (R) is written in Morse oscillator 
form, and hence analytical Morse oscillator states are 
used for X„(-R), while the Xe,l(R) are integrated numer- 
ically and normalized such that (see Child, Appendix A 



/ X E >,L(R)XE,L(R)dR = S(E-E>). 
Jo 



(A6) 
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